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The spin dynamics in [TbPc2]'^ and [DyPc2]'' single molecule magnets have been investigated 
by means of muon and nuclear spin-lattice relaxation rate measurements. The correlation time 
for the spin fluctuations was found to be close to 0.1 ms already at 50 K, about two orders of 
magnitude larger than the one previously found in other lanthanide based single molecule magnets. 
In [TbPc2]°two different regimes for the spin fluctuations have been evidenced: a high temperature 
activated one involving spin fluctuations across a barrier A ~ 880K separating the ground and 
first excited states and a low temperature regime involving quantum fluctuations within the twofold 
degenerate ground-state. In [DyPc2]''a high temperature activated spin dynamics is also evidenced 
which, however, cannot be explained in terms of a single spin-phonon coupling constant. 
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Molecular nanomagnets have attracted major atten- 
tion in the last decade owing to their possible techno- 
logical applicability as logic units in computer*^ or in 
life sciences, as contrast agents^. So far, the most stud- 
ied molecular nanomagnets are the ones based on tran- 
sition metal ions, which are characterized by a twofold 
degenerate ground-state separated from the first excited 
states by a barrier induced by the magnetic anisotropy^. 
The stability of the ground-state and the applicability of 
these compounds above liquid helium temperature (T) is 
nevertheless prevented by the small magnetic anisotropy 
which characterize transition metal complexes^. In or- 
der to increase the anisotropy the synthesis of lanthanide 
based molecular magnets has been envisaged and re- 
cently the double-decker phthalocyaninato lanthanide 
(Ln) complexes [Pc2Ln]^ (Pc= C32H16N8) have been 
recognized as potential single molecule magnets (SMM) 
which could function above cryogenics temperatures^'^. 
In fact, [Pc2Tb]~ and [Pc2Dy]~ complexes are the first 
mononuclear lanthanide systems exhibiting a very slow 
relaxation and quantum tunneling of magnetization^i^i^. 
It has been recently observed^ by means of nuclear 
magnetic resonance (NMR) that in [TbPc2]"TBA+ x 
N [TEA] Br, where N represents the degree of dilution 
within the organic tetrabutylammonia (TBA) matrix, 
the characteristic correlation time (tc) for the spin fluctu- 
ations becomes of the order of a fj,s at liquid nitrogen T,^ 
much longer than in any other molecular magnet. These 
relatively long Tc originates from the huge barrier sep- 
arating the twofold degenerate \ J — 6,m — ±6) crystal 
field (CP) ground-state from the first excited | J = 6, to = 
±5) state. In [TbPc2]"TBA+ x N [TBA]Br compounds, 
barriers as high as A ~ 920 K have been found^. The 
CP splitting is determined both by [TbPc2]~ structure 
and by the arrangement of the other molecules around 
it. In fact, upon varying N different values for A have 
been estimated. Moreover, it has been observed that at 
temperatures T <C A quantum fluctuations, related to 
the tunneling between to = ±6 levels, arise. 

In order to increase the stability of the ground-state the 



CP splitting has to be enhanced and the tunneling rate 
lowered. In this work we have considered the effect of the 
one-electron oxidation of [Pc2Ln] ~ . The resulting neutral 
[Pc2Ln]'^ complex has an unpaired tt electron which res- 
onates between the Pc molecules, yielding a compression 
of the cage around Ln"^"*" ion and possibly to an increase 
in the CP splittin g^^i^^i^'^ . Hereafter we present an in- 
vestigation of the spin dynamics in two different neutral 
compounds, [TbPc2]''and [DyPc2]°, by means of NMR 
and muon spin relaxation (/xSR). The T dependence of 
Tc has been derived by means of muon and nuclear spin- 
lattice relaxation rates and in [TbPc2]" found to be close 
to 0.1 ms, already at 50 K, much longer than the one 
found in [TbPc2]-TBA+ x N [TBA]Br^. The T depen- 
dence of the relaxation rate is determined by the spin 
fluctuations induced by spin-phonon coupling. Accord- 
ingly, information on the CP level splitting and on the 
spin-phonon coupling was derived. 

[TbPc2]''and [DyPc2]'' complexes were synthesized ac- 
cording to the templatingiS reaction of the phthaloni- 
trile precursor in boiling amyl alcohol in presence of 
DBU under inert atmosphere for 24 h. Treatment of 
[Ln(acac)3xn H2O] with phthalonitrile gave both forms, 
the anionic [Pc2Ln]^ and the neutral [Pc2Ln]", as ob- 
served in other earlier reports^^. The crude reaction mix- 
ture was presorbed on active (H2O-0%) basic alumina 
oxide. This treatment helps to convert unstabilized an- 
ionic form to the neutral one. Purification was carried 
out by column chromatography on basic alumina oxide 
(deactivated with 4.6% H2O, level IV) with chloroform 
- methanol mixture (10:1) as eluent. Analytically pure 
powder samples were achieved by additional radial chro- 
matography on silica gel followed by recristallization from 
chloroform - hexane mixture. 

NMR measurements have been performed by us- 
ing standard RP pulse sequences. While above 100 K 
the full NMR line could be irradiated and the spectra 
obtained from the Pourier transform of half of the echo, 
below 40 K the spectra became rather broad (Pig. la) 
and they could be derived from the envelope of the echo 
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FIG. 1: a) NMR spectrum in [TbPca]" at T = 19 K (cir- 
cles), compared to the much narrower spectrum obtained at 
T = 202 K (line). The linewidth at half intensity decreases 
from a few MHz at 19 K to about 37 kHz at 202 K. b) Time 
evolution of the muon polarization in [TbPc2]'' sample nor- 
malized to its value for i ^ at sbc selected temperatures. 



amplitude upon varying the irradiation frequency. This 
broadening is a clear evidence of the slowing down of 
the spin dynamics in these SMM. Further information 
in this respect can be obtained from the study of the 
nuclear spin-lattice relaxation rate 1/Ti. 1/Ti was 
estimated from the recovery of the nuclear magnetiza- 
tion after a saturating RF pulse sequence. The recovery 
was found to be a stretched exponential, which evidences 
a distribution of relaxation rates for the nuclei, likely 
to be related to differences in the hyperfine coupling be- 
tween Ln'^+ and the 16 protons of each Pc molecule. The 
T dependence of 1/Ti, reported in the inset of Fig. 2, 
evidences a progressive increase upon cooling from RT 
down to T ~ 140 K, while for T < 40 K 1/Ti is found 
to progressively decrease. In the intermediate T range 
140 K> r > 40 K the observation of NMR signal is 
prevented by the short relaxation times. 

In order to explore the spin dynamics in these interme- 
diate T range we have studied [TbPc2]"and [DyPc2]°by 
means of /iSR, which allows to measure much faster re- 
laxation rates and to work at low fields, thus perturbing 



the SMM only to a minor extent. Zero (ZF) and longi- 
tudinal field (LF) /iSR experiments were carried out at 
ISIS pulsed muon facility on MUSR beam line. The back- 
ground asymmetry contribution was found to be around 
10% out of the 30% total initial asymmetry, for both 
samples. The time decay of the muon polarization Pf^it) 
in the samples shows a different behaviour for tempera- 
tures above and below T* ~ 90 K, for [TbPc2]", and ~ 60 
K, for [DyPc2]°(Fig. lb). Above T* the decay follows a 
stretched exponential behaviour Pfj,{t) — Aexp{—{Xt)^), 
with an exponent (3 ~ 0.5 and A ~ 20%. 
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FIG. 2: T dependence of the muon longitudinal relaxation 
rate in [TbPc2]°and [DyPc2]°for H = 1000 Gauss (circles) 
and for H — 6000 Gauss (squares). The lines are the best fits 
according to Eq. 3. In the inset the T dependence of ^H l/Ti 
is reported for the two compounds, for H = 1 Tesla. 

Below T* a marked decrease of the initial asymmetry is 
found. This decrease has to be associated with the onset 
of very low-frequency fluctuations. In particular, if the 
hyperfine field at the muon fluctuates with a character- 
istic correlation time Tc 3> 1/7/iV < A/i^ >, with 7^ the 
muon gyromagnetic ratio and V < A/i^ > the root mean 
squared field distribution probed by the muons, the decay 
of Pfi{t) is given by the static Kubo-Toyabe functioni^, 
which in zero-field is characterized by a fast initial decay 
and by the subsequent recovery of the muon polarization 
to A/3. At ISIS the initial fast decay cannot be detected, 
owing to the pulsed nature of the muon source and, ac- 
cordingly, a loss in the initial polarization is observed 
and only the long-tim e tail detected. At longitudinal 
fields H ~ V< A/i2 > an increase in the tail amplitude 
is expected. Since a LF of 1 kGauss (Fig. lb) is found to 
cause an increase of Pp(t ~^ 0)/P^(G) to a{H) ~ 0.55 for 
[TbPc2]", one can estimate V< A/i^ > ~ 650 Gauss^. 
Thus, one has to expect that belo w T* Ln'^+ moments 
are frozen over a timescale (l/j^V < A/i^ >) ~ 10 ns. 
Still, even within this frozen configuration low-energy 
excitations are present and drive muon spin-lattice re- 
laxation. Hence, for T < 80 K, one has that Pp(i) = 
a{H)Aexp{-{\t)f^). 



3 



The T dependence of A, derived by fitting the decay 
of the muon polarization according to the procedure re- 
ported above is reported in Fig. 2. One notices a peak 
around 90 K for [TbPc2]°and around 60 K for [DyPca]", 
in agreement with the NMR findings. The intensity of 
the peak is observed to scale with the inverse of the field 
intensity, a situation found when the frequency of the 
fluctuations is close to Larmor frequency lol- 
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FIG. 3: Muon and spin-lattice relaxation rates are re- 
ported as a function of inverse T in a semi- logarithmic plot, 
in order to evidence the activated behaviour for T > T* . The 
lines are the best fits yielding A = 881 K for [TbPc2]°and 
A = 610 K for [DyPc2]°. 

From the T dependence of the spin-lattice relaxation 
rates it is possible to obtain information on the charac- 
teristic correlation time for the spin fluctuations, on the 
energy barrier between the CF ground-state and the first 
excited states and on spin-phonon coupling. In view of 
the energy difference between the muon (nuclear) hyper- 
fine levels and the m levels of Ln'^"'" , direct muon (nuclear) 
relaxation processes involving an electron spin excitation 
are forbidden. The effective relaxation processes are in- 
direct ones, involving a muon (nuclear) spin fiip without 
change in |m|. This is possible as the dipolar hyper- 
fine Hamiltonian contains terms coupling the transverse 
components of the hyperfine field hx,y to Jj. Owing to 
spin-phonon scattering processes each CF level is charac- 
terized by a finite life-time (r^) which yields a Lorentzian 
broadening and A (or 1/Ti) can be written aai^ 
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with (A/i'^) the mean-square amplitude of the hyperfine 
field fluctuations, Em the eigenvalues of the CF levels and 
Z the corresponding partition function. The life-time for 
the m levels can also be expressed in terms of the tran- 
sition probabilities Pm,m±i between m and m ± 1 in the 
form I/tto = Pm,rn+i + Pm,m-i- Siuce the life-times are 
mainly determined by spin-phonon scattering processes 



those can be expressed in terms of the CF eigenvaluesi^ 

{Em±l — Em)^ 



,/3(i5„±i-£;„) „ I 



(2) 



with C spin spin-phonon coupling constant and (3 = 1/T. 
Since in [TbPc2]° A = £^^=±5 - Em=±6 > ksT over aU 
the explored T range, Eq. 2 can be simplified in the form 
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with Tc — exp{A/T)/{CA^). From the amplitude of A at 
T* one can estimate a root mean square amplitude for 
the fiuctuating field at the muon in [TbPc2]°-\/ (A/i-^) ~ 
770 Gauss. This value is close to the one of the field 
distribution giving rise to the fast initial relaxation of the 
muon polarization, described by Kubo-Toyabe function. 
By fitting the data with Eq. 3 it is possible to estimate a 
spin-phonon coupling constant C ~ 3000 Hz/K'^, of the 
same order of magnitude of that found in other SMM^ii^. 
It is noticed that for T > 100 K, TcOJl <C 1 and from Eq. 
3 one has A c>c exp{A/T). In fact, by plotting cither l/Ti 
or A data vs. 1/T in a semi-logarithmic scale (Fig. 3), 
one finds a nice linear behaviour, consistent with a value 
of A ~ 880 K. 

On the basis of Eq. 3, in the light of the data reported 
in Fig. 2, it is possible to derive the T dependence of Tc 
for [TbPc2]°. In [TbPc2]° Tc shows two regimes (Fig. 4): 
a high T activated one and a low T one, for T < 50 K, 
where the correlation time is constant. The high T be- 
haviour describes spin fluctuations among m — ±6 and 
m = ±5 levels, while the low T trend rather signals tun- 
neling processes among m = +6 and m = —6 levels. The 
total fluctuation rate is given by the sum of the fluctua- 
tion rates associated with each process and accordingly 
one has 
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(-), = CA3e-^/^ + (l), 

Tc Tc 



(4) 



with a tunneling rate (I/tc)* ~ 90/is. It is noticed (Fig. 
4) that the T-dependence of Tc can be reproduced very 
well with the above equation. It is now useful to compare 
NMR and fiSR results in [TbPc2]°with AC susceptibility 
ones^^ii^ii^. In fact, also these measurements evidence 
a high T activated behaviour with energy barriers which 
appear to be only slightly smaller than the ones found 
here. On the other hand, the role of tunneling processes 
in the AC susceptibility measurements has not been ad- 
dressed for the Pc2Ln complexes and just the frequency 
dependence of the AC susceptibility peak discussediiii^. 
Nevertheless, we remark that also the imaginary part of 
the AC susceptibility is characterized by a low T flatten- 
ing which might originate from tunneling processesii*i^. 
Still, in comparing these techniques one should consider 
that local techniques as NMR and /iSR are sensitive to 
the spin fluctuations of each molecule even if they do not 
yield a net variation in the macroscopic magnetization, 
namely are silent in the AC susceptibility—. 
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FIG. 4: T dependence of the correlation time for the spin 
fluctuations in [TbPc2]°and [DyPc2]°(for T > T*) , derived 
from A data reported in Fig. 2, on the basis of Eq. 3. The 
solid lines are the best fits according to Eq. 4. 

Now, if one considers the T dependence of A in 
[DyPc2]°(Fig. 2), one notices a behaviour quite similar 
to the one found in [TbPc2]'^, with a maximum around 
60 K characterized by a sUghtly lower amplitude. This 
decrease in the amplitude has to be associated with a 
lower hyperfine coupling, corresponding to a root mean 
square amplitude for the fluctuating field at the muon 
y/ (A/i^) ~ 590 Gauss. Since Dy'^^ magnetic moment is 
slightly larger than Tb^+ one, this decrease rather indi- 
cates a slightly larger average distance between the rare 
earth ion and the muon in [DyPc2]'^. In view of the simi- 
larity with [ThPcal^one could try to fit also [DyPc2]°A(r) 
data with Eq. 3. One finds a good fit only for T > 40 K 



and one can estimate a barrier for the spin fiuctuations 
around 610 K (Fig. 3) and a spin-phonon coupling con- 
stant C ~ 1780 Hz/K^. However, one notices that for 
T < 40 K the data are not so well reproduced. In fact, 
for [DyPc2]*'it is likely that the two lowest doublets are 
separated by an energy barrier in the tens of K ranged 
and that the high barrier derived from A(T) corresponds 
to the one between the first and second excited doublets 
|J, ±m). Hence, in order to appropriately fit the data 
one should resort to Eq. 1. However, even by allowing 
different possible values for the CF splittings a poor fit 
of [DyPc2]°data is obtained. It is likely that one should 
consider that the spin-phonon coupling constant is not 
the same for the scattering processes yielding high and 
low energy transitions. 

In conclusion, by means of fiSK and NMR spin-lattice 
relaxation rates we have derived the T dependence of the 
characteristic correlation time for the spin fiuctuations 
in [TbPc2]°and [DyPc2]°SMM. In [TbPc2]°, although the 
barrier between the ground and first excited states is close 
to the one of [TbPc2]" SMM diluted in a TEA matrix^, 
Tc is found to be almost two orders of magnitude larger 
for T < 50 K. This suggests that in the neutral SMM the 
structural arrangement around the SMM itselfiS, leads to 
an enhancement of the CF symmetry and to a decrease 
of the quantum fiuctuations within the twofold degener- 
ate ground-state. Thus, it appears of major relevance, in 
order to enhance Tc and to be able to employ SMM at 
economically affordable temperatures, to carefully inves- 
tigate the effect of the modifications in the SMM coordi- 
nation on the spin dynamics. 

Useful Discussions with F. Borsa are gratefully ac- 
knowledged. The research activity in Pavia was sup- 
ported by Fondazione Cariplo 2008-2229 research funds. 
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